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ABSTRACT 


This  document  presents  a  general  explicit  solution 
to  the  powered  flight  dynamics  of  a  rocket  vehicle 
having  constant  thrust  and  constant  effective  exhaust 
velocity,  and  demonstrates  the  use  of  this  solution 
in  guidance  techniques. 
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1.  INTRODU  •  .  ION 


In  the  latter  part  of  1959,  the  author  developed  an  explicit  solution  to  the 
powered  flight  dynamics  of  a  rocket  vehicle  having  constant  thrust  and  constant 
effective  exhaust  velocity.  This  solution  was  uaed  as  the  basis  for  the  guidance 
equations  for  the  Me rcury/ Atlas  flights.  Since  that  time,  guidance  equations 
based  on  this  solution  have  been  developed  for  use  on  Project  Gemini  and 
Ranger/Mariner,  as  well  as  for  other  programs  of  a  more  classified  (security) 
nature . 

This  explicit  solution  and  its  application  to  guid;.nce  equations  was  documented 
in  Reference  (1).  ‘  Since  that  time,  many  people  have  examined  this  document 
as  well  as  other  techniques  for  guiding  a  continuously  burning  vehicle  to  pre¬ 
scribed  values  of  altitude  and  velocity  at  burnout.  As  a  result  of  these  investi¬ 
gations  \t  appears  clear  that  for  the  class  of  vehicles  under  consideration: 

(1)  any  guidance  technique  which  accurately  controls  burnout  altitude 
and  velocity  vector  must  be  based  on  an  explicit  solution  of  the 
powered  flight  dynamics. 

(2)  the  solution  which  was  originally  developed  and  documented  in 
Reference  (1)  is  essentially  definitive  (the  details  of  mechanization 
will  of  course  depend  to  a  cenain  extent  on  particular  mission 
requirements ). 

Although  the  notation  in  this  document  differs  from  that  used  in  Reference  (1) 
there  are  no  non-trivial  changes  in  the  subject  matter.  Some  additional  infor¬ 
mation  or  explanatory  material  has  been  provided  for  various  peripheral  a»  as 
because  those  omissions  have  resulted  in  considerable  confusion  by  many 


"""An  Explicit  Method  of  Guiding  a  Vehicle  from  an  Arbitrary  Initial  Position  and 
Velocity  to  a  Prescribed  Orbit,”  by  D.  MaePherson,  dated  13  February  1961, 
Report  No.  TPR -594(1 565-0 1  )TN- 1 . 
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readers  of  Reference  (1).  The  order  of  presem  .on  in  the  Basic  Vehicle 
Kinematics  section  has  been  changed  for  greater  clarity  (the  manner  of  presen¬ 
tation  used  here  was  suggested  by  Mr.  M.  Tangora).  Various  small  errors  in 
Reference  (1)  were  discovered  by  a  number  of  people;  these  have  hopefully  been 
corrected.  Mr.  W.  Brocato  has  been  especially  helpful  in  this  regard. 

This  explicit  solution  to  the  powered  flight  dynamics  of  a  rocket  vehicle  is  of 
fundamental  importance  for  guidance  equations  on  sophisticated  missions; 
because  the  total  effect  of  control  commands  is  to  determine  the  position  and 
velocity  of  the  vehicle  at  that  point  where  ability  to  control  terminates  (it  is 
obvious  that  only  certain  restricted  values  of  position  and  velocity  can  be 
obtained).  The  values  of  position  and  velocity  that  are  desired  or  required  at 
control  termination  (thrust  termination)  may  be  specified  in  a  variety  of  ways. 

On  an  ICBM  mission,  for  example,  it  is  possible  to  make  the  ve.'ocity  at  thrust 
termination  a  function  of  position  at  thrust  termination,  and  thereby  eliminate 
the  need  for  position  control.  On  the  other  hand,  burnout  position  and  velocity 
are  independently  specified  on  many  space  missions.  The  explicit  solution 
contained  lierein  can  be  utilized  to  attain  these  desired  conditions  at  thrust 
termination  regardless  of  the  origin  of  the  specification  (the  specification  may 
be  considered  another  part  of  the  guidr-.nce  equations). 

In  the  development  it  is  assumed  that  vehicle  thrust  attitude  changes  (or  control) 
take  place  in  two  mutually  perpendicular  planes  which  arc  conventionally  desig¬ 
nated  as  pitch  and  yaw.  It  is  assumed  that  the  thrust  magnitude  cannot  be  con¬ 
trolled  (except  by  termination). 

This  document  has  been  written  in  as  general  a  manner  as  possible  with  emphasis 
placed  on  the  principles  involved  rather  than  on  the  sophistications  and  details 
of  mechanization  which  might  be  desirable  when  using  this  solution  in  a  detailed 
set  of  guidance  equations.  For  example,  the  time  difference  between  generated 
control  commands  and  effective  time  of  data  and/ or  effective  time  of  command 
enactment  should  be  accounted  for;  but  this  has  not  been  discussed  because  a 
bookkeeping  technique  of  this  type  is  not  conceptually  difficult  or  fundamental  to 
the  system. 
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The  equivalent  relationship 


T1 


1  - 


c* 


(tj  -  t) 


(2.4) 


will  be  more  convenient  at  times. 

The  following  dimensionless  variable  is  of  fundamental  importance. 


U  =  - 


-  t) 


(2.5) 


Note  that 


(2.6) 


f 

The  constant  value  of  U  is  in  a  sense  a  measure  of  the  vehicle's  dynamic  charac¬ 
teristics.  Then  "time  to  go"  until  t  =  tj  can  be  found  from 


(2.7) 


The  evaluation  of  Ty  will  be  discussed  after  the  immediately  following  develop¬ 
ment  of  an  appropriate  expression  for  use  in  the  evaluation  of  U. 

The  acceleration  when  U  =  0  is  a^j  end  will  be  redefined  as  ay  Then  from 
Equations  (2.3),  (2.5),  and  (2.6) 


ay  =  Cc’t' 


_  _ 

®T  -rrTT 


(2.8) 


(2.9) 
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> 

The  general  equation  of  motion  along  V  .will  now  be  developed. 
The  general  vehicle  force  diagram  is  shown  in  Figure  2.  1. 


Fy  *  mo  y 


Figure  2.  1 


These  forces  will  be  divided  into  two  components;  one  parallel  to,  and  one  per¬ 
pendicular  to  velocity.  The  latter  component  is  of  no  interest  in  the  immediate 
discussion  since  it  ciianges  only  the  direction  and  not  the  magnitude  of  the 
velocity.  The  force  balance  along  V  may  be  written  as: 


cos  (P  ••  r)  -  g  sin  r  + 


cos  ^ 


Defining 


1 

cos  (P  -  n  +  g  sin  r  -  -pjj-  cos  9 


(2. 10) 
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Then 
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Then 


I 


The  magnitude  of  V'  is  still  unknown.  However,  the  derivative  at  any  time  is 

Lj 

simply 


or  numerically 


If  V'  were  known  initially,  the  above  equation  could  be  used  to  find  proper 

Lj 

succeeding  values.  This  technique  will  be  used  as  defined  by  the  following 
equations; 


Vf  -  V  1 


U  ^  1  -  c 


c>:< 


whe  re 


'’h  =  '^L  ,  -  “l-' 

n  n- 1 


a,  is  from  Equation  (2.  10) 

Lj 

T  is  length  of  computation  cycle 


Vy  is  a  preselected  constant 


(2.13) 


(2.14) 
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The  value  of  may  be  arbitrarily  choHcn;  however, 

reasonably  r'*alistio  value  so  that  Vj^  is  near  7.  -o  vhen 
from  Equation  (2.  13)  that  when  V  = 


it  is  desirable  to  use  a 
V  =  Vj.  It  is  obvious 


V 


Uj  =  1  -  e 


Lf 

c* 


and  that  elapsed  time  until  V  =  is 


U,.  -U 

“IT" 


The  uncertainty  in  is  reflected  in  U^.  in  the  above  mechanization.  Although 

the  above  mechanization  is  valid,  the  following  alternative  has  generally  been 
used  in  the  past. 

I  he  time  difference  6Ty  biitween  the  points  where  V  =  Vj  and  V  =  is 


6T 


U 


(2.15) 


where  5^  is  the  average  thrust  acceleration  between  times  tj  and  tj.  Using 
Equation  (2.9) 


A]  +  6Ty 

/  a,,,  dt 

A,  1 


6T 


U 


=  a 


in  (I  -  U6Ty) 


U 


■U6T 


'2. 16) 


U 


1  -  U 


Then  the  time  interval  until  velocity  equals  deotred  velocity  is  (see  Figure  2.2) 


T 


+  6Tu 


(2. 17) 


U  »  0 


Figure  2.  2 


The  time  interval  Ty  can  be  mechanized  as  (from  Equation  2.7) 


Since  U  can  be  found  by  differencing  successive  values  of  U.  Ty  can  be  deter 
mined  from  knowledge  of  c'!'  and  measurements  of  velocity  only.  (See  Appendix  A 
for  effect  of  error  in  knowledge  of  c*.  )  It  should  especially  be  noted  that  while 
the  thrust  must  be  constant,  it  need  no-  be  known.  Since  U  is  linear  in  time 
(from  Equation  (2.5)),  it  is  an  appropri.  le  (quantity  on  which  to  apply  smoothing, 
and  smoothing  should  be  applied  on  U  in  preference  to  an  equivalent  nonlinear 
quantity  (V  for  example). 
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This  intogral  must  be  predicted. 


It  is  obvious  that  the  magnitude  of  will  not  be  known  exactly  when  T^,  is 
large;  but  as  becomes  smaller  the  value  of  the  integrand  in  Equation  (2.  18) 
not  only  becomes  more  predictable  but  also  requires  prediction  for  shorter 
periods  of  time.  Techniques  which  have  been  used  in  this  prediction  are  dis¬ 
cussed  in  Appendix  D. 


The  characteristics  of  Ty  and  6Ty  which  have  been  cited  lead  to  the  conclusion 
that  the  percentage  error  in  T^,  is  always  small.  Simulation  has  verified  this 
conclusion  by  consistently  demonstrating  that  errors  in  Tj^,  are  less  than  (and 
usually  considerably  less  than)  one  per  cent  in  the  absence  of  errors  in  position 
and  velocity  (see  simulation  results  in  Section  4).  It  should  be  noted  that  the 
separation  of  the  thrust  effects  and  the  "loss"  effects  (into  and  or  the 
equivalent)  is  necessary  in  a  practical  mechanization  in  order  to  achieve 
stability.  This  separation  provides  a  decoupling  which  causes  first  order 
perturbations  in  loss  estimates  to  have  only  second  order  effects  on  time  to 
go  (Tj,)  estimates.  If  decoupling  is  not  provided,  perturbations  in  loss  esti¬ 
mates  corrupt  the  estimate  of  thrust  acceleration  which  causes  large  variations 
in  time  to  go.  A  combination  of  small  loss  perturbations  and  heavy  smoothing 
on  time  to  go  can  prevent  instability,  but  this  is  not  a  desirable  mechanization 
for  obvious  reasons. 


Note  that  Equation  (2.  9)  gives  thrust  acceleration  in  terms  of  the  same  smoothed 
quantities  that  determine  T^j,  and  that  this  thrust  acceleration  can  be  extrapolated 
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backwards  or  forwards  in  time  with  no  more  error  (percentagewise)  than  exists 
in  the  current  estimate.  In  particular,  the  acceleration  at  engine  cutoff  is 


‘U 


1  -  U6T 


U 


(2. 19) 


If  thrust  and  structural  mass  of  the  vehicle  are  known,  Equation  (2.  19)  can  be 
used  to  predict  the  mass  of  propellant  remaining  at  engine  cutoff. 

Figure  ?..  3  is  a  collection  of  equations  from  this  section  which  are  needed  for 
computation  of  time -to -go  and  steering  commands.  It  is  not  suggested  that  this 
is  the  exact  set  of  equations  which  would  be  mechanized.  Note  that  since  U  is 
used  in  the  solution  as  well  as  U,  the  j.  .  ohanlcatlon  of  Eq.  {?,.  13)  should  be  a 
truncation  of  the  exponential  series  instead  of  a  curve  fit. 

It  will  be  convenient  to  note  for  future  reference  that  the  total  integral  of  thrust 
acceleration  between  times  t  and  t^  is 


AV 


.■f' 


=  V, . 


v  +  V,  -  V 


Lf 


(2.  20) 
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Figure  2.  3 
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3.  EFFECTS  OF  E  ATTITUDE 


Equations  relating  missile  attitude  profile  and  vehicle  position  and  velocity  will 
be  developed  in  this  section.  These  relationships  can  then  be  used  for  control 
purposes.  Pitch  plane  dynamics  will  be  discussed  first  as  they  are  the  most 
general.  Vaw  dynamics  will  then  be  a  special  case  of  this  general  solution. 

Great  mathematical  simplicity  would  result  if  the  development  could  be  made  in 
an  inertial  cartesian  coordinate  system.  The  approach  taken  in  the  following 
derivation  will  be  to  solve  the  problem  in  such  an  inertial  frame  and  to  add 
terms  to  the  solution  to  account  for  the  differences  between  this  model  and 
reality.  In  order  to  validly  use  this  approach  it  is  necessary  to  eliminate  the 
anomalies  caused  by  the  following  items  in  the  physical  environment: 

(1)  The  "true"  inertial  coordinate  system  is  at  the  geocentroid  and  not 
rotating,  while  the  desired  vehicle  position  is  most  conveniently 
specified  as  a  radius  from  the  geocentroid. 

(2)  The  existence  of  a  gravitational  field. 

Item  (1)  can  be  compensated  for  by  use  of  the  following  coordinate  frame  (see 
Figure  2.1): 

Let  V  =  inertial  velocity  of  the  vehicle 

r  =  distance  from  the  geocentroid  to  the  vehicle 

Then 


V 

sin  r  = 


(3.  1) 
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Since  these  variables  are  derived  relative  to  r,  ley  "rotate"  with  an  anpular 

V  cos  r  ’ 

rate  =  - - -  but  they  are  truly  inertial  quantities  since  they  are 

measured  in  an  inertial  coordinate  system.  (The  rotation  is  "stopped"  while 

the  computation  is  being  made.  ) 

Likewise  a  pitching  rate 


-  V  cos  r 
r 


(3.2) 


is  necessary  to  keep  thrust  attitude  (p^)  constant  (relative  to  r).  This  pitching 
rate  will  be  transmitted  in  addition  to  any  pitching  rate  defined  for  other 
purposes. 


The  gravitational  and  coriolis  accelerations  in  these  coordinates  produce  a  com¬ 
bined  acceleration 


.  [  V  cos  r]^ 

^Ep  ~  ^r  r 


(3.3) 


where  g  is  the  component  of  gravitational  attraction  along  r  and  may  be  defined 

^  _  > 
in  any  desired  degree  of  complexity  (the  component  perpendicular  to  r  is  of  no 

consequence  here,  and  is  negligible  in  any  case).  The  effect  of  g^p  can  be 

compensated  for  by  defining  a  missile  attitude 


(3.4) 


and  adding  this  attitude  to  the  attitude  desired  for  position  and  velocity  control. 
This  specifically  assumes  that 


sin  P  =  P 


(3.5) 


-14- 


a  rnlationship  that  will  be  assumed  for  si*i  city  throughout  the  development. 
While  this  approximation  is  valid  for  most  practical  trajectories,  less  approxi¬ 
mate  relationships  are  easily  derived  (see  Appendix  H).  A  pitching  rate 


b> 

gP 


(3.6) 


must  be  transmitted  to  maintain  the  correct  value  of 
Equation  (3.4)  is  given  in  Appendix  C. 


The  differentiation  of 


The  equations  governing  position  and  velocity  in  an  inertial  plane  will  now  be 
developed.  These  vehicle  attitudes  must  satisfy  the  following  relationships, 
using  Equation  (3.5): 


(3.7) 


(3.8) 


A  necessary  condition  to  the  solution  of  Equations  (3.  /)  and  (3.8)  is  that 

P  =  A  +  A,f,(t,  -  t) 

^mp  o  1  1'  f  ' 

where  A^  and  A^  are  arbitrary  constants  and  fj  is  an  arbitrary  function  of 
(tj  -  t),  and,  of  course 

“mp  "  ^  ['l“f  ■  *’] 
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The  functional  form  of  f  -  t)  is  restricted  on*  .1  that  it  must  satisfy 
Equations  (3.7)  and  (3.8);  and,  in  theory,  need  be  neither  analytic  nor  continuous. 
There  are,  however,  several  rather  obvious  physical  reasons  for  restricting  the 
physical  form  of  f  j(lj  “  t)*  The  attitude  must  be  continuous  and  attitude  rates 
should  be  small  on  practical  trajectories.  These  requirements  are  satisfied  by 
setting 


fj(tj  -  t)  t^  -  t 


or 


“  *0  ■*  *i<‘t 


t) 


This  makes  w  constant,  with  the  resulting  mathematical  advantages, 
rnp 

Several  developments  have  been  advanced  by  different  writers  with  various 
physical  assumptions  on  the  problem  of  maximizing  performance;  the  results  of 
all  argue  that  cither  0  or  sin  0  or  tan  0  should  be  linear  in  time.  The 
fact  that  these  "proofs"  are  not  strictly  applicable  to  the  real  case  is  often  over- 
look<;tI  or  not  emphasized,  but  they  do  indicate  that  a  constant  attitude  rate  solu¬ 
tion  should  be  relatively  efficient  (particularly  when  g^-p  is  small  relative  to  a,j.). 
This  constant  attitude  rate  solution  is  generally  used  and  will  be  developed  here. 
Other  types  of  solutions  can  be  easily  developed  along  the  same  lines.  A  develop¬ 
ment  employing  two  constant  attitudes  has  been  used  on  Ranger/Mariner. 

It  is  possible  to  choose  a  constant  attitude  0^  which  would  satisfy  Equation  (3.7) 
[without,  in  general,  satisfying  Equation  (3.8)].  For  this  attitude  Equation  (3.'  / 
reduces  to 


V 

po 


+  ?  AV 
P  e 


(3.9) 
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s 


Then  missile  attitude  ran  be;  expressed  as 


ft  =  ^  +  P  I-  -  '*»  (t,  -  t) 

mp  p  prf  pr  f 


(3.  11) 


Equation  (3.  10)  satisfies  Equation  (3.7)  by  requiring 


I’prf  -  -  ‘’j"'  '  “ 


(3.  12) 


The  notation  w  is  adopted  to  signify  that  it  is  needed  only  to  provide  position 
P  ^ 

(radius)  control.  The  implications  and  advantages  of  this  formulation  of  Pj^p 
will  be  discussed  following  the  development.  The  various  P  quantities  are  shown 
schematically  in  Figure  3.  1. 


Figure  3.  1 


Expressions  for  Pp^f  and  Wp^  will  be  found  by  su  .tituting  Equation  (3.  11)  into 
Equation  (3.8)  and  combining  with  Equation  (3.  12).  This  yields  (see  AppendixB) 


u 

pr 


■ 

E _ E_ 


(°6  - 


'°6^E 


pr 


w  lie  re 


(3.  13) 


(3.  14) 


(3.  15) 


The  numerator  in  Equation  (3.  13)  is  the  altitude  "error"  that  would  exist  if  Wp^ 

were  held  at  zero;  the  denominator  is  the  appropriate  scaling  factor  between 

0^  and  this  altitude  error, 

pr 

The  desired  present  missile  attitude  is  the  sum  of  the  various  attitude  compo¬ 
nents  which  have  been  generated  above,  so 


w  T_  -  p 
pr  E  pg 


-  (1 


-  C,)w  T„ 
6  pr  E 


(3.  16) 


The  actual  missile  attitude  pp  which  is  determined  by  any  suitable  method  will 
not  in  general  be  equal  to  P^  (although  in  "steady  state"  the  difference  should  be 
quite  small)  so  that  there  will  be  an  additional  transient  pitch  turning  rate 


V  ■  ‘‘’o  -  Pp^^ 


(3.  17) 
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where  G  is  a  gain  function.  Then  the  total  pitch  turning  rate  is 


(j)  =0)  +  CJV  +  U  +0) 

P  Pg  P  Pr  pt 


(3.  18) 


After  an  initial  trahsient,  w  will  be  very  small;  and  w  and  w*  are  always 
'  Pf  Pg  P 

small  because  of  the  physics  of  the  problem.  The  magnitude  of  is  in  effect 

a  measure  of  the  feasibility  of  satisfying  the  specified  burnout  conditions,  since 
if  an  excessive  amount  of  steering  is  required  it  may  not  be  practical  or  even 
possible  to  satisfy  these  burnout  conditions  (this  is  why  Equation  (3.  5)  is  usually 
satisfactory).  Note  that  limiting  the  magnitude  of  will  modify  the  burnout 
altitude  requirement  without  affecting  the  ability  to  satisfy  the  velocity  require¬ 
ment  in  any  way  whatsoever.  This  feature  is  especially  convenient  near  cutoff, 
since  there  is  practically  no  capability  to  change  altitude  and  the  inclusion  of  a 

limit  on  w  allows  the  equations  to  function  in  a  rational  manner  without 
pr 

necessitating  any  other  changes. 

It  should  be  noted  for  completeness  that  the  superposition  which  has  been  assumed 
is  not  strictly  valid  for  non-zero  values  of  V^.  This  effect  is  extremely  small 
for  practical  thrust  levels  (burning  times),  vanishes  as  T^,  goes  to  zero,  and 
will  not  be  discussed  for  these  reasons. 

Figure  3.2  is  a  collection  of  equations  from  this  section  which  are  needed  for 
pitch  steering.  As  in  Figure  2.3,  these  equatio.us  are  not  necessarily  ideal  for 
mechanization. 

Yaw  steering  may  be  mechanized  in  a  number  of  ways.  It  would  be  possible,  for 
example,  to  define  a  desired  orbit  plane  with  displacement  y  and  a  yaw  velocity 
Vy,  perpendicular  to  this  plane.  Substitution  of  appropriate  notation  in 
Equations  (3.  7)  through  (3.  18)  is  all  that  is  required  to  make  these  equations 
applicable  to  yaw  steering  for  this  case.  Equations  (3.2),  (3.4)  and  (3.6)  are 
replaced  by 


(3  =  w  =  w>l‘  =  0 

yg  yg  y 


a  condition  that  exists  for  any  form  of  yaw  steering. 
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(3.19) 


It  may  be  unnecessary  and  thus  undcsii  ’  xe  (from  a  peiformance  standpoint)  to 
control  the  orbit  plane  with  the  precision  available  in  the  above  equations.  A 
somewhat  simpler  technique  would  be  to  choose  a  point  in  the  desired  orbit  plane 
(a  "pseudo-target")  and  require  that  the  intersection  of  the  actual  orbit  plane 
and  the  desired  orbit  plane  contain  this  point.  This  can  be  insured  by  defining  a 
yaw  velocity  (V^)  perpendicular  to  the  plane  containing  the  missile  (r),  the 
target  (r.^),  and  the  geocentroid;  and  then  requiring  that  be  zero  at  cutoff. 

This  may  be  mechanized  by  defining 


V 

y 


>  > 
r.j.  X  r  . 

"r^r  sin 


> 

V 


(3.20) 


where  is  the  angle  between  r.j.  and  r.  Requiring  only  that  =  0  is  in  effect 
setting  the  terni  corresponding  to  equal  to  zero,  and  thus  replacing 
Equation  (3.  16),  (3.  17)  and  (3.  18)  by 


u' 


(Py  -  Py)G 


(3.21) 


At  cutoff  the  vehicle  will  be  displaced  from  the  desired  burnout  plane  by 
approximately 


^  cutoff 

Ay  =  /  V  dt  (3.22) 

•'liftoff  y 

The  integrand  in  Equation  (3.22)  vanishes  at  liftoff  and  cutoff,  and,  in  tl  , 
absence  of  malfunction,  never  becomes  very  large.  Since  the  vehicle  is  in  the 
desired  plane  at  the  pseudo-target,  this  pseudo-target  can  for  many  missions  be 
selected  such  that  the  magnitude  of  Ay  is  of  little  or  no  importance.  If  the  mag- 
nif'ide  of  Ay  is  sufficiently  small  in  the  presence  of  normal  dispersions  (as  it 
usually  is)  the  simplicity  and  fuel  economy  resulting  from  this  pseuoo-target 
technique  justify  its  use  for  these  missions. 
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4.  SIMULATIOI  .ESULTS 


Guidance  equations  utilizing  the  preceding  explicit  solution  have  been  simu¬ 
lated  by  several  people  for  use  in  achieving  a  variety  of  burnout  conditions  (both 
with  and  without  noise)  for  a  variety  of  missions.  In  the  absence  of  errors  in 
positions  and  velocity  information,  approximate  errors  of  injection  into  an  orbit 
(if  c‘l=  is  known)  are  given  by  the  following  tabulation: 


Parameter 
Velocity  Magnitude 
Velocity  Orientation 
Altitude 


3ir  dispersion 
.  1  foot  per  second 
.5(10”^)  radian  =  .  003  degree 
150  feet 


Injection  at  a  relatively  large  P  does,  for  rather  obvious  reasons,  increase  the 
altitude  dispersions,  although  this  dispersion  still  remains  small,  These  dis¬ 
persions  are  due  to  data  lag,  response  lag,  smoothing  lag  (filters  are  used  even 
in  the  absence  of  noise),  variability  of  thrust  and  c*  from  assumed  constant 
value  due  to  the  use  of  influence  coefficients  in  the  simulation  of  engine  thrust 
and  mass  flow,  etc.  In  an  actual  mission  noise  and  various  other  hardware  type 
effects  (such  as  cutoff  impulse  uncertainty)  will  cause  considerably  larger  injec¬ 
tion  errors.  Some  of  these  errors  can  be  reduced  by  the  use  of  a  vernier. 

Lack  of  knowledge  of  the  exact  value  of  will  introduce  a  bias  in  time-to-go 
only  if  smoothing  is  introduced  (see  Appendix  I).  The  size  of  this  bias  will 
depend  not  only  on  the  "error"  in  c*  but  also  on  the  "time  constant"  of  the 
smoothing  filter;  however,  reasonable  values  of  these  quantities  produce  a  ve'*'' 
small  effect  on  velocity  magnitude  and  negligible  effect  on  altitude  and  velocity 
orientation.  Errors  due  to  noise  will  not  be  discussed  except  to  note  that  the 
extreme  precision  (linearity)  of  the  "time-to-go"  solution  minimizes  the 
severity  of  problems  arising  from  this  area  (absolute  linearity  would  give 
optimum  results). 
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APPENDIX  A 

EFFECTS  OF  ER^.  JRS  IN  c« 


Differcntiatinp  Equation  (A.  2) 


Then 


Equation  (A.  7)  id  tabulated  below: 


m 

0 

my 

C*!' 

1 

0 

1.5 

0. 199 

2 

0.307 

3 

0.451 

4 

0.538 

5 

0.598 

10 

0.744 

oo 

1 

(A.  7) 


It  is  apparent  that  9Ty/9c’i'  is  satisfactorily  small  at  all  times,  and  that  the 
percentage  error  vanishes  as  goes  to  zero. 
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APPEND1>.  B 


INTEGRAL  EVALUATION 


The  notation  will  be  changed  in  the  mathematical  development  of  this  appendix 
to  provide  greater  simplicity  and  clarity.  Final  results  will  be  transformed 
into  the  original  notation.  The  changes  are  as  follows: 

1.  The  subscript  o  is  used  to  denote  values  "now"  and  unsubscripted 
values  are  used  to  denote  values  between  "now"  and  final. 

2.  The  time  axis  is  defined  so  that  t  =0. 

o 

From  Equation  (2.4)  the  final  acceleration  is 


1 


a  c’!' 


o 


(B.  1) 


so  that 


c*  _  c* 

aftf  a^tj 


The  following  relationships  will  be  useful 


(E.2) 


(B.3) 
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Then 


which  in  combination  with  Equatio;  {B.2)  and  (B.4);  and  with  the  definition 


AV 
_ c 

a^t^ 


(B.5) 


becomes 


=  c-tj.05 


(B.6) 


From  Equation  (B.  1) 


a 

o 
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which  in  combination  with  Equations  (B.?‘  and  (B,5) 


gives 


-  05] 


(B.7) 


Also 


•'O  Jo  ■ 


,t  dt  dt  -c* 


dt 


•  i r/o^r  o*  o*]  =  [- ^“5]  =  [- °5  *  s(^)] 


Defining 


O,  ? 


c«0^ 


6  ■  AV 


(B.8) 


Then 


X  “*  *  =  '■*'/  [-  I  +  °5  "  06''  -  °5>] 


(B.9) 


Integrating  Equation  (3.  12) 


dt  -  0 
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Combining  with  Equations  (D.4)  and  {B.7) 


The  numerical  values  of  and  T^,  are  id  nt-cal,  and  making  this  substitution 
in  Equations  (B.  11)  and  (B.  10)  produces  Equations  (3.  13)  and  (3.  14)  respectively. 
Also 


Q6  = 


c*  c*!* 
^  "  "f'^E 


(B.  12) 
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APPENDIX  C 


From  Eqiuitiori  (2.3) 


Then 


Neglecting  small  terms 

d(3 

_ S 

dt 

or  with  Equation  (3.2) 


^Ed 

w  =  — -£  'Vj  2u>« - 

pg  dt  p 


(C.4) 


It  may  be  easily  shown  that 


dt  1  ^c*:'  J 


(C.5) 


The  fact  that  di*j  / dt  is  extremely  small  is  fortunate  for  obvious  reasons. 

a 
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Then 


+  r  (t  -  t^)l  dt 


/ '  *  V ‘  '  'X '  '^1 


or 


r 


(^) 


(D.3) 


Defining 


X  5  1  .  cos  (P  -  r)  % 


ID.4) 


Then  with  Equation  (D.  2) 


Pf 


J  'a^^dt  =J^^^a^Xdt=|jr^^a^[(Pj-rj)^  +  ZCP^-F^HP'  -r)(t-tj)  +  (p'  -  T' )^(t  -  tj)^]dt 

zj\^\dt  =  (Pj  -  r^)^^^a^dt-2{pj-rjKP’  -r')jr^^a^tj-t)dt  +  (p'  -r')^  ^a^(tf-t)^dt 
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Using  conventions  of  Appendix  B 


2  jT  a^\dt  =  |(pf  -  -  2(Pf  -  r^)(P'  -  r')t^  +  (p'  -  r)^t/]  j  a^dt  + 

+  [2(P^  -  rj)(p'  -  r)  -  (p*  -  r')^2tjj  + 

+  (P'  -  r')2  I  a  t^dt 

•'o  ^ 

From  Eqs.  (B.  1),  (B.  4),  (B.  7),  (B.  8),  (D.  2)and  (D.4) 
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With  Eq.  (3.  4) 


J 

/ 


f  ,2 

dt  .^^2 

- - =  C*tj 

I  -  -St 


From  Eq.  (B.  2).  (B.  5)  and  (B.  7) 


0*5*  ,  c* 

=  1  + 


'Cf 


a  t, 
o  f. 


=  1  + 


c«(l  -  Q3) 
AVi 


c*  .  /  c*' 


=  (i  .Q^r +  2(1  -03)^+  - 


Then 


^  ,2 

/a  t 

o  ‘ 


dt  =  c’i't 


AV  y  P  jjl  C  ^ 

^(l  .Q^)^4  2(l  .Q^)+^  — 

c  c 


-  AV  tf 
e  f 


Defining 


Q.^  s  (P-  -  r')(tj  - 1)  =  (Pj  -  Fj)  -  (p  -  r) 


1 

2 


(D.6) 


(D.  7) 
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Then  with  Eq.  (D.6)  and  (0.7)  Eq.  (D.5)  I 


.ornes 


\dt  AV 


X  +  (p  -  r)  (I  -  o^)  + 


Then 


j 

/ 


a  dt  =  XAV 
LX  e 


(D.S) 


Whe  re 


x  =  x+(p  -r)o^(i  -  o^)  +  -j"- 


<>  -  °6>'  r 


c* 


Experience  shows  that  the  following  approximation  gives  very  accurate  results 


X  = 


(D.9) 


where  is  given  in  Equation  (D.7). 
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AJtiiough  aoroclynamic  forms  could  be  handled  in  a  similar  manner, 
tori('s  normally  considered  it  is  very  satisfactory  to  set 


'  “i-A  =  " 


Then  Equation  (2,  18)  can  be  written  as 


^Lf  = 


where  \  is  defined  in  Equation  (D.  9). 


(^) 


AV  \ 
e 


for  trajec- 

(D.  10) 


(D.  11) 
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APPENP’X  E 


EVALUATION  OF 


From  Equation  (2.  16) 

06T  (U6T..)^  (U6T„)^ 

l+— ^  + + .  (E.l) 

From  Equation  (2.9) 


in(l  -  06Ty) 


'U 


■U6T 


U 


'1  1 

"  1  -  U6Tu 


+  UoTy  f  (IJ6Ty)^  +  (U6Ty)^  +  .  .  . 


ao  that 


Then 


ay  +aj 
~  2  ■  “ 


tr6T,,  (06T„)^  {06T,.)^ 

1  +  — __ii  +  - _bL_  +  - -iL_  + 

Zee 


‘U 


(06Ty)^ 
- g - 


+  - 


(U6Ty)- 


(E.2) 


It  is  obvious  irom  Equation  (E.2)  that  the  approximation 


(E.3) 
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is  very  pooc!.  It  can  be  easily  shown  that  (for  c*  10,000  ft/sec)  the  error  in 
injection  velocity  incurred  by  using  Equation  {r  ,  is  approximately 


where  6V  and  are  in  units  of  ft/ sec.  Usually  «  800  ft/ sec. 
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APPENDIX  F 


EFFECTS  OF  COMPUTATION  CYCLE  LENGTH 
AND  AUTOPILOT  RESPONSE 


The  injection  errors  in  velocity  magnitude  and  yaw  velocity  are  very  small  (less 
than  1  ft/sec.  )  for  all  computing  cycle  lengths  less  than  10  sec.  Pitch  velocity 
errors  as  a  function  of  major  computation  cycle  time  have  also  been  determined 
from  simulation  and  are  shown  in  Figure  F-1.  Position  errors  are  correspond¬ 
ingly  small.  These  error  numbers  do  not  include  any  errors  except  those 
arising  from  the  guidance  equations  themselves.  The  hardware  contribution  to 
the  injection  errors  will  be  relatively  independent  of  computation  cycle  but  will 
be  somewhat  larger  for  the  longer  computing  intervals.  It  is  assumed  in  all 
cases  that  the  determination  of  position  and  velocity  information  will  be  done 
frequently  enough  to  prevent  computation  error  build-up.  For  the  longer  compu 
tation  cycle  times,  this  may  require  the  up-dating  of  position  and  velocity  infor¬ 
mation  in  sub-cycles.  It  should  be  noted  that  the  errors  shown  in  Figure  F-1 
are  to  a  certain  extent  due  to  lack  of  sophistication  in  the  mechanization  of  w* 

and  w  (these  quantities  are  not  constants  during  the  computation  cycle).  These 
1  ^  Is 

errors  could  be  greatly  reduced  by  the  use  of  offsets  or  a  more  sophisticated 
mechanization. 


Figure  F-1 
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1  he  effects  of  lack  of  knowledge  of  autopilot  '•esoonses  have  been  approximately 
simulated  on  the  IBM  7090.  The  simulation  w..ich  was  used  has  a  unity  autopilot 
with  a  delay  of  any  desired  magnittide.  The  guidance  equations  assume  a  delay 
of  Z  seconds.  Simulations  were  run  with  actual  delays  on  the  dynamics  varying 
from  U  to  6  sec.  which  is,  of  course,  far  outside  any  possible  lack  of  knowledge 
of  response  characteristics. 

When  errors  in  pitch  velocity  occur  due  to  very  long  computing  cycles  or 
because  autopilot  response  is  grossly  misjudged,  pitch  steering  commands  can 
assume  large  or  even  limited  values  near  vehicle  thrust  termination.  These 
turning  commands  as  well  as  the  other  steering  commands  used  during  the  flight 
are  not  oscillatory  and  in  this  sense  are  not  unstable.  The  existence  of  this 
stability  has  been  not  only  proven  by  simulation,  but  may  be  deduced  as  follows 
from  past  experience  and  knowledge  of  system  behavior. 

In  a  radio  guidance  system,  vehicle  thrust  attitude  must  be  deduced  from  radar 
position  and  velocity  information.  The  radar  noise  normally  produces  estimates 
of  attitude  that  are  quite  ina.ccurate,  but  these  attitude  estimates  are  then 
smoothed  in  a  digital  filter  with  a  sufficiently  long  time  constant  to  permit  a 
reasonably  accurate  attitude  to  be  obtained.  This  technique  works  very  well  (as 
has  been  demonstrated  on  Project  Mercury)  even  when  considerable  noise  is 
present  in  the  radar  data.  When  an  inertial  system  is  used,  there  may 
corresponding  "noise”  in  a  position  and  velocity  mcasur' d  by  the  IMU.  These 
errors  result  from  rotational  dynamics  and  from  the  fact  that  the  IMU  is  not 
Iccated  at  the  center  of  gravity  of  the  vehicle  as  v/ell  as  quantization  of  measure¬ 
ments.  Although  these  errors  are  considerably  smaller  than  the  corresponding 
noise  errors  in  a  radio  guidance  system,  the  technique  for  handling  thorn  can  be 
the  same;  i.e.  ,  raw  attitudes  are  o  btained  from  position  and  velocity  information 
and  smoothed  appropriately.  Since  the  "noise"  errors  on  an  inertial  system  are 
much  smaller  than  those  for  a  radio  guidance  system,  the  problems  arising  from 
this  area  are  correspondingly  smaller.  Since  the  only  data  required  by  the 
guidance  equations  are  position  and  velocity  information,  the  above  technique 
will  insure  stability. 
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In  some  guidance  techniques  n.>t  based  on  .  n  xplicit  powered  flight  solution  it 
is  not  pos'  ible  or  practical  to  smooth  thrust  attitude,  and  this  can  lead  to  insta¬ 
bility  under  certain  circumsta  ices  when  these  less  sophisticated  techniques  are 
used. 
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APPEN^'Ty  H 


LESS  APPROXIMATE  RELATIONSHIPS 


More  general  equations  which  pa  tially  remove  the  approximation  of  Equation 
(3.  5)  will  be  given  without  discussion. 

Equation  (3.4)  becomes 


sin  P 


gP 


(H.  1) 


Equation  (3,  10)  becomes 


(H,2) 


Equation  (3.  16)  becomes 


sin  P  =  sin  P  + 
o  P 


sin  P  -  (1  - 
gP 


Q/-)w  T_ 
6'  pr  E 


(H.3) 


Equation  (3.  17)  becomes 


w  .  -  I  sin  P  -  sin  P  ]G 
pt  o  p^ 


(H.4) 


Equation  (3.  18)  becomes 


W  =  («: 

P  F 


[u>  -f  a'  +  w 
Pg  pr  pt 

cos  P 


(H.5) 
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Eqviation  (G.2)  becomes 


P*  =  sin 
P 


V  -  V  , 
pn  pn-1 

aj-r 


+  P 


gP 


-  Q 


2n 


(H.6) 
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NOTAi  '■  J 


The  notation  has  been  subdivided  into  symbols  and  subscripts.  Some  of  the 
symbols  appear  only  in  combination  with  subscripts. 


SYMBOLS 


g 


m 

r 


t 

T, 


"U 

V 

u 


r 

T 


-  Thrust  acceleration 

-  Average  thrust  acceleration  over  the  time  interval  ^iT^j 

-  Effective  exhaust  velocity 

-  Aerodynamic  force 

-  Thrust 

-  The  total  force  exerted  on  the  vehicle  by  gravitational  fields 
divided  by  the  mass  of  the  vehicle 

-  Mass  of  the  vehicle  including  unexpended  propellants 
•  Distance  from  the  goocentroid  to  the  vehicle 

-  Time 

-  Time  to  go  until  the  desired  velocity  is  attained 

-  Time  to  go  until.  U  “  0 


-  Inertial  velocity  of  the  vehicle 

-  Fundamental  dimensionless  variable  -  defined  in  Equation  (2.5), 
computed  from  Equation  (2.13) 

-  Missile  attitude 

-  Component  of  thrust  attitude  necessary  for  position  and  velocity 
control  in  an  inertial  cartesian  coordinate  frame 

-  Flight  path  angle  (positive  up  from  local  horizontal) 

-  Length  of  computational  cycle 
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SYMBOLS  (continued) 

-  Vehicle  attitude  turning  rate 

-  See  Equation  (3.2) 

SUBSCRIPTS 

f  -  Final  (when  velocity  equals  desired  velocity) 

g  -  Associated  with  gravity 

p  -  Pitch 

r  -  Associated  with  geocentric  radius 

t  -  Transient 

y  -  Y  aw 

E  -  Effective 

L  -  Associated  with  acceleration  "loss'* 
o  -  Piesont  value 

'i  '  Thrust 

U  -  Related  to  U 

?.  -  Associated  with  angle  of  attack 
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